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Abstract 

Using measured radial velocity data of nine double lined spectroscopic binary systems 
NSV 223, AB And, V2082 Cyg, HS Her, V918 Her, BV Dra, BW Dra, V2357 Oph, and YZ 
Cas, we find corresponding orbital and spectroscopic elements via the method introduced by 
Karami & Mohebi (2007a) and Karami & Teimoorinia (2007). Our numerical results are in 
good agreement with those obtained by others using more traditional methods. 
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1 Introduction 

Determining the orbital elements of binary stars helps us to obtain fundamental information, 
such as the masses and radii of individual stars, that has an important role in understanding 
the present state and evolution of many interesting stellar objects. Analysis of both light and 
radial velocity (hereafter RV) curves, derived from photometric and spectroscopic observations, 
respectively, yields a complete set of basic absolute parameters. One historically well-known 
method to analyze the RV curve is that of Lehmann-Filhes (cf. Smart, 1990). In the present 
paper we use the method introduced by Karami & Mohebi (2007a) (=KM2007a) and Karami & 
Teimoorinia (2007) (=KT2007), to obtain orbital parameters of the nine double- lined spectro- 
scopic binary systems: NSV 223, AB And, V2082 Cyg, HS Her, V918 Her, BV Dra, BW Dra, 
V2357 Oph, and YZ Cas. 

The NSV 223 is a contact system of the A type and the mass ratio is believed to be small 
(Rucinski et al. 2003a, b). The large semiamplitudes ,Ki, suggest that the orbital inclination 
is close to 90°. The spectral type is F7V and the period is 0.366128 days (Rucinski et al. 
2003a, b). The AB And is a contact binary. The spectral type is G8V. The period is 0.3318919 
days. It is suggested that observed period variability may be a result of the orbital motion 
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in a wide triple system. The third body should then have to be a white dwarf (Pych et al. 
2003,2004). V2082 Cyg is most probably an A- type contact binary with a period of 0.714084 
days. The spectral type is F2V (Pych et al. 2003,2004). In HS Her, the effective temperatures 
were found to be Ti = 15200 it 750K and T2 = 7600 it 400-ftr for the primary and secondary 
stars, respectively (Cakirli et al. 2007). The secondary component appears to rotate more 
slowly. The presence of a third body physically bound to the eclipsing pair has been suggested 
by many investigators. The two component are located near the zero-age main sequence, with 
age of about 32 Myr (Cakirli et al. 2007). It is classified as an Algol- type eclipsing binary 
and single-lined spectroscopic binary. The spectral type of more massive primary component 
is B4.5V. The effective temperature is about 15200 it 750K for the primary component and 
7600 lb 400i^r for the secondary component. The period is 1.6374316 days (Cakirli et al. 2007). 
The V918 Her is an A-type contact binary. The spectral classification is A7V. The period of 
this binary star is 0.57481 days (Pych et al. 2003,2004). The BV Dra and BW Dra have a 
circular orbit. Prom spectrophotometry the components of BV Dar are classified as F9 and F8 
while the components of BW Dra are G3 and GO. The period of BV Dra is 0.350066 days and 
for BW Dra is 0.292166 days (Batten & Wenxian 1986). The V2357 Oph was classified as a 
pulsating star with a period of 0.208 days. The spectral type is G5V (Rucinski et al. 2003a,b). 
The spectral type of primary component of YZ Gas is AlV and for the secondary is F7V. The 
period of this binary is 4.4672235 days (Lacy 1981). 

This paper is organized as follows. In Sect. 2, we give a brief review of the method of 
KM2007a and KT2007. In Sect. 3, the numerical results are reported, while the conclusions are 
given in Section 4. 



2 A brief review on the method of KM2007a and KT2007 

One may show that the radial acceleration scaled by the period is obtained as 

-2ttK / ,,Z 

PZ = 



-/VTA . / _^ A X 

2^sm(^cos (_-ecosa;)j 



X 1 1 -I- e cos ^ — a; -I- COS ^(^ — ecosa;)^| , (1) 

where the dot denotes the time derivative, e is the eccentricity, oj is the longitude of periastron 
and Z is the radial velocity of system with respect to the center of mass. Also 

2tt asini 

where P is the period of motion, a is the semimajor axis of the orbit and inclination i is the 
angle between the line of sight and the normal of the orbital plane. 

Equation (1) describes a nonlinear relation, PZ = PZ{Z,K,e,uj), in terms of the orbital 
elements K, e and w. Using the nonlinear regression of Eq. (1), one can estimate the parameters 
K, e and u, simultaneously. Also one may show that the adopted spectroscopic elements, i.e. 
mp/nis, nipSm^ i and rUs sin^i, are related to the orbital parameters. 



3 Numerical Results 

Here we use the method of KM2007a and KT2007 to derive both the orbital and combined 
elements for the nine different double lined spectroscopic systems NSV 223, AB And, V2082 
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Cyg, HS Her, V918 Her, BV Dra, BW Dra, V2357 Oph, and YZ Cas. Using measured radial 
velocity data of the two components of these systems obtained by Rucinski et al. (2003a,b) for 
NSV 223 and V2357 Oph, Pych et al. (2003,2004) for AB And, V2082 Cyg, and V918 Her, 
Cakirli et al. (2007) for HS Her, Batten & Wenxian (1986) for BV Dra and BW Dra, Lacy 
(1981) for YZ Cas, the fitted velocity curves are plotted in terms of the photometric phase in 
Figs. 1 to 9. 

Figures 10 to 27 show the radial acceleration scaled by the period versus the radial velocity 
for the primary and secondary components of NSV 223, AB And, V2082 Cyg, HS Her, V918 
Her, BV Dra, BW Dra, V2357 Oph, and YZ Cas, respectively. The solid closed curves are 
results obtained from the non-linear regression of Eq. (1), which their good coincidence with 
the measured data yields to derive the optimized parameters K, e and uj. The apparent closeness 
of these curves to ellipses is due to the low, or zero, eccentricities of the corresponding orbits. 
For a well-defined eccentricity, the acceleration- velocity curve shows a noticeable deviation from 
a regular ellipse (see Karami & Mohebi 2007b). 

The orbital parameters, K, e and ui, obtaining from the non-linear least squares of Eq. (1) 
and the combined spectroscopic elements including nip sin^ i, nis sin'^ i, {ap + a^) sin i and mp/mg 
obtaining from the estimated parameters K, e and ui for the nine systems are tabulated in Tables 
1 to 18, respectively. The velocity of the center of mass, V^m, is obtained by calculating the 
areas above and below of the radial velocity curve. Where these areas become equal to each 
other, the velocity of center of mass is obtained. Tables 1 to 18 show that the results are in 
good accordance with the those obtained by Rucinski et al. (2003a,b) for NSV 223 and V2357 
Oph, Pych et al. (2003,2004) for AB And, V2082 Cyg, and V918 Her, Cakirh et al. (2007) for 
HS Her, Batten & Wenxian (1986) for BV Dra and BW Dra, and Lacy (1981) for YZ Cas. 

4 Conclusions 

Using the method introduced by KM2007a and KT2007, we obtain both the orbital elements 
and the combined spectroscopic parameters of nine double lined spectroscopic binary systems. 
Our results are in good agreement with the those obtained by others using more traditional 
methods. In a subsequent paper we intend to study the other different systems. 
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Figure 1: Radial velocities of the primary and secondary components of NSV 223 plotted 
against the photometric phase. The observational data have been deduced from Rucinski et 
al. (2003a,b). 
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Figure 2: Same as Fig. 1, for AB And. The observational data have been measured by Pych et 
al. (2003,2004). 
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Figure 3: Same as Fig. 1, for V2082 Cyg. The observational data have been derived from Pych 
et al. (2003,2004). 
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V918 Her 
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Figure 5: Same as Fig. 1, for V918 Her. The observational data have been derived from Pych 
et al. (2003,2004). 
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Figure 6: Same as Fig. 1, for BV Dra. The observational data have been derived from Batten 
& Wenxian (1986). 
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Figure 7: Same as Fig. 1, for BW Dra. The observational data have been derived from Batten 
& Wenxian (1986). 
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Figure 8: Same as Fig. 1, for V2357 Oph. The observational data have been derived from 
Rucinski et al. (2003a,b). 
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Figure 9: Same as Fig. 1, for YZ Cas. The observational data have been derived from Lacy 
(1981). 
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NSV 223 (primary) 




Figure 10: The radial acceleration scaled by the period versus the radial velocity of the primary 
component of NSV 223. The solid curve is obtained from the non-linear regression of Eq. (1). 
The plus points are the experimental data. 
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Figure 11: Same as Fig. 10, for the secondary component of NSV 223. 
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Figure 12: Same as Fig. 10, for the primary component of AB And. 
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Figure 13: Same as Fig. 10, for the secondary component of AB And. 
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Figure 14: Same as Fig. 10, for the primary component of V2082 Cyg. 
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Figure 15: Same as Fig. 10, for the secondary component of V2082 Cyg. 
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Figure 16: Same as Fig. 10, for the primary component of HS Her. 
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Figure 17: Same as Fig. 10, for the secondary component of HS Her. 
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Figure 18: Same as Fig. 10, for the primary component of V918 Her. 




13 



BV Dra (primary) 




Figure 20: Same as Fig. 10, for the primary component of BV Dra. 
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Figure 21: Same as Fig. 10, for the secondary component of BV Dra. 
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Figure 22: Same as Fig. 10, for the primary component of BW Dra. 
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Figure 23: Same as Fig. 10, for the secondary component of BW Dra. 
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Figure 24: Same as Fig. 10, for the primary component of V2357 Oph. 
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Figure 25: Same as Fig. 10, for the secondary component of V2357 Oph. 
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Table 1: Orbital parameters of NSV 223. 



This Paper Rucinski et al. (2003a,b) 



Primary 






Vcm {kms~'^) 


-21.31 ±0.98 


-21.66(2.15) 




42.56 ±0.03 


40.39(2.64) 




0.004 ± 0.001 




o 


279.40 ± 12.04 




Secondary 






Vcm {kms~^) 


-21.31 ±0.98 


-21.66(2.15) 


Kg {kms~^) 


300.76 ± 0.04 


297.98(4.18) 


e 


e.s = Rp 











Table 2: Spectroscopic elements of NSV 223. 

Parameter This Paper Rucinski et al. (2003a,b) 

nip sin^i/MQ 1.3447 ± 0.0007 — 

m,sSm^i/MQ 0.1903 ± 0.0002 — 

{ap + as) sin i/Rq 2.4834 ± 0.0005 — 

ms/mp 0.1415 ± 0.0001 0.136(10) 

{mp + ms)sm^i/MQ 1.535 ± 0.001 1.47(9) 



Table 3: Same as Table 1, for AB And 





This Paper 


Pych et al. (2003,2004) 


Primary 






Vcm (kms^^) 


-27.26 ± 0.66 


-27.53(0.67) 


Kp (kms^^) 


232.69 ± 0.02 


232.88(0.83) 




0.00109 ± 0.00005 






230.69 ± 3.22 




Secondary 






Vcm {kms~^) 


-27.26 ± 0.66 


-27.53(0.67) 


Kg {kms~^) 


132.43 ± 0.01 


130.32(1.17) 




Cg — Gp 




In 


iOs = uOp- 180° 





Table 4: Same as Table 2, for AB And 



Parameter This Paper Pych et al. (2003,2004) 

mp sin^ i/Mq 0.6071 ± 0.0001 — 

rus sin^ i/Mq 1.0668 ± 0.0002 — 

{ap + as)smi/RQ 2.3943 ± 0.0002 — 

mp/ms 0.5691 ± 0.0001 0.560(7) 

{mp + ms)sin^i/MQ 1.6739 ± 0.0004 1.648(20) 
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Table 5: 


: Same as Table 1, for V2082 Cyg 


This Paper Pych et al. (2003,2004) 


Primary 






Vcm {kms~^) 


-34.19 ±0.39 


-34.12(0.58) 


Kp (kms~^) 


32.254 ± 0.003 


33.16(0.51) 


e 1 


0.0008 ± 0.0005 






266.48 ± 0.06 










Vcm {kms~^) 


-34.19 ±0.39 


-34.12(0.58) 


Kg (kms~^) 


139.50 ± 0.01 


139.38(0.99) 


e 


65 = €p 






Us = iOp- 180° 




Table 6: 


: Same as Table 2, for V2082 Cyg 


Parameter 


This Paper 


Pych et al. (2003,2004) 


mpSin^i/MQ 


0.3045 ± 0.0001 




rus s\v? i/MQ 


0.07039 ± 0.00002 




{ttp + Us) sini/i?0 


2.4232 ± 0.0002 




ms/mp 


0.23121 ± 0.00004 


0.238(5) 


{mp + nis) sin^ z/Mq 


0.3749 ± 0.0001 


0.380(7) 


Table 7: Same as Table 1, 


for HS Her 


This paper Cakirli et al. (2007) 


Primary 






Vcm {kms~^) 


-14.24 ±0.87 


-12.8 ± 1.9 


Kp [kms~^) 


91.15 ±0.02 


93.4 ±2.8 




0.048 ± 0.001 


0.05 ± 0.01 




268.66 ± 0.03 




Secondary 






Vcm {kms~^) 


-14.24 ± 0.87 


-12.8 ± 1.9 


Kg (kms"^) 


232.13 ± 0.01 


239.4 ± 3.2 




(^s — (^p 


0.05 ± 0.01 


In 


iOs = iOp- 180° 




Table 


8: Same as Table 2, for HS Her 


Parameter 


This study 


Cakirh et al. (2007) 


rup sin^ I/Mq 


4.116 ±0.001 




nis sin^ i/MQ 


1.616 ±0.001 




{ttp + as) sini/R 


10.458 ± 0.001 


10.76 ± 0.50 


ms/mp 


0.3927 ± 0.0001 


0.39 ± 0.05 
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Table 9: Same as Table 1, for V918 Her 





This Paper 


Fycn et ai. (zOUo,zU04j 


Primary 






Vcm {kms~^) 


-25.76 ± 0.78 


-25.72(0.74) 


Kp {kms~^) 


199.54 ± 0.01 


199.37(1.72) 




0.0002 ± 0.0001 




In 


315.14 ± 12.17 




Secondary 






Vcm {kms~^) 


-25.76 ± 0.78 


-25.72(0.74) 


Kg {kms~^) 


53.577 ± 0.001 


53.93(0.49) 




6s = 6p 




In 


ujs = ujp- 180° 





Table 10: Same as Table 2, for V918 Her 



Parameter 


This Paper 


Pych et al. (2003,2004) 


mpSiiL^i/MQ 


0.20442 ± 0.00002 




rus s\v? i/Mq 


0.7613 ±0.0001 




mp/rris 


0.26851 ± 0.00001 


0.271(5) 


{nip + iTis) s\v? I/Mq 


0.9657 ± 0.0001 


0.968(21) 



Table 11: Same as Table 1, for BV Dra 





This Paper 


Batten & Wenxian (1986) 


Primary 






Vcm (kms^^) 


-61.23 ±0.77 


-65.2 ± 1.2 


Kp {kms~^) 


94.119 ±0.002 


93.9 ± 1.4 




0.0025 ± 0.0003 




In 


256.06 ± 11.18 




Secondary 






Vcm {kms~^) 


-61.23 ±0.77 


-57.7 ± 1.7 


Kg (kms"^) 


139.50 ± 0.01 


139.38(0.99) 




65 = Cp 




In 


Us=ojp- 180° 





Table 12: Same as Table 2, for BV Dra 



Parameter This Paper Batten &; Wenxian (1986) 



apSini/m^ 


0.453 ± 0.001 


0.452 ± 0.007 


agSini/lO^ 
nip sin^ i/MQ 


1.1248 ±0.0001 


1.124 ±0.010 


0.9106 ± 0.0002 


0.911 ± .020 


nig sin^ i/Mq 


0.36679 ± 0.00006 


0.366 ±0.010 


nig /nip 


0.40278 ± 0.00004 


0.402 ± 0.007 
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Table 13: Same as Table 1, for BW Dra 





A- liiO J. a^JKDL 


BflttpTi Rt- WpnviflTi nQSR"! 


Primary 






Vcm [kms 


-61.14 it 0.73 


-64.2 ± 0.9 


Kp [kms ) 


TA cr\ _L n no 
7(J.59 ± (J. (Jo 


70. ± 1.1 




Bp = Bg 




In 


270.58 ± 0.44 




Secondary 






Vcm {kms~^) 


-61.14 ±0.73 


-57.4 ± 1.4 


Kg {kms~^) 


252.018 ± 0.005 


251.5 ± 1.7 




0.0026 ± 0.0001 




In 


u;s = Up- 180° 




Table 14: Same as Table 2, for BW Dra 


Parameter 


This Paper 


Batten & Wenxian (1986) 


Up sinz/10^ 


0.2836 ± 0.0001 


0.283 ± 0.004 


as sini/10^ 


1.01249 ± 0.00002 


1.010 ±0.007 


mpSin^ i/Mq 


0.7939 ± 0.0001 


0.791 ± 0.015 


rus s\v? i/Mq 


0.2224 ± 0.0001 


0.222 ±0.005 


mp/iris 


0.2801 ± 0.0001 


0.280 ± 0.005 


Table 15: Orbital parameters of V2357 Oph. 




This Paper 


Rucinski et al. (2003a,b) 


Primary 






Vcm {kins~^) 


-17.79 ± 0.63 


-19.12(1.64) 


Kp [kms"^) 


44.29 ± 0.07 


44.12(1.63) 




Cp — Cg 




In 


227.12 ±0.23 




Secondary 






Vcm {kms~^) 


-17.79 ±0.63 


-19.12(1.64) 


Kg {kms~^) 


192.52 ±0.02 


190.93(2.90) 




0.0134 ± 0.0004 




In 


LOg = u;p- 180° 





Table 16: Spectroscopic elements of V2357 Oph. 



Paraiii(;t;C}r 


Tliis Pajjcr 


Rudnt,ki vX al. (2003a,b) 


mpSiiL^i/MQ 


0.4647 ± 0.0004 




rUg sin^ i /Mq 


0.1069 ± 0.0002 




mg/rup 


0.2300 ± 0.0004 


0.231(10) 


{rup + mg) s\v? i/Mq 


0.5716 ± 0.0007 


0.560(32) 
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Table 17: Same as Table 1, for YZ Cas 





This Paper 


Lacy (1981) 


Primary 






Vcm {kms~^) 


7.78 ± 0.58 


8.14 ±0.06 




74.039 ± 0.002 


73.35 ± 0.21 




0.0034 ± 0.0001 


0.0 ± 0.003 


In 


270.078 it 0.395 




Secondary 






Vcm {kms~^) 


7.78 ± 0.58 


8.14 ±0.06 


Kg {kms~^) 


124.49 dr 0.05 


125.7 ±0.5 




— Gp 




In 


ujs = ijOp- 180° 





Table 18: Same as Table 2, for YZ Cas 
Parameter This Paper Lacy (1981) 



CpSini/lO^ 
Qg sini/lO^ 

rUp/MQ 



4.548 ± 0.003 
7.6470 ± 0.0001 
2.274 ± 0.002 
1.352 ±0.001 



4.508 ± 0.012 
7.727 ± 0.031 
2.31 ±0.01 
1.35 ±0.02 
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